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Abstract. Gaussian Klauder coherent states are discussed in the context of the infinite well 
quantum model, otherwise known as the particle in a box. A supersymmetric partner system 
is also presented, as well as a construction of coherent states in this new system. We show 
that these states can be chosen, in both systems to have many properties usually expected for 
coherent states. In particular, they yield highly localised wave packets for a short period of 
time, which evolve in a quasi-classical manner and which saturate approximately Heisenberg 
uncertainty relation. These studies are elaborated in one- and two-dimensional contexts. Finally, 
some relations are established between the gaussian states being mostly used here and the 
generalised coherent states, which are more standardly found in the literature. 


1. Introduction 

The study of coherent states in quantum mechanics is a subject very well documented 
[H El El 01151E!- These states have been approached from different points of view and we can 
cite, for example, their definition as eigenstates of an annihilation operator and as minimium 
uncertainty states. They have also been defined for a large amount of solvable quantum systems 
and their properties have been analysed 

One approach, that has not deserved much attention until now, is the one using the so-called 
gaussian Klauder coherent states j3]. They are given as a superposition of energy eigenstates 
that leads to a good localisation in the phase space of the system under consideration. For the 
harmonic oscillator, these last states are a good approximation of the usual coherent states. It 
has been shown that it is also the case for the infinite well and the Morse systems 00 main]. 
It means that they are good states for investigating physical properties of the systems under 
consideration (behaviour of the position and momentum observables and uncertainty relation, 
for example). 

Even if most of the works were focussing on one dimensional (ID) systems, some results have 
been given in the two dimensional (2D) case |10| . In particular, the 2D infinite well, having 
a quadratic degenerate energy spectrum, has attracted some attention and the construction of 
coherent states had to be adjusted noma . 

All these considerations have been extended using a supersymmetric (SUSY) approach [13]. 
Indeed, SUSY partner Hamiltonians have been constructed in ID and 2D using intertwining 
relations and factorisation methods na na ng. These new Hamiltonians are shown to be 


deeply related to the original ones. The coherent states systems are as well closely related to 
the original ones. 

In this work, we consider the infinite well quantum system both in ID and 2D and their SUSY 
partners. We then construct the Gaussian Klauder coherent states (GCS) and investigate their 
properties. In Section 2, we show how the SUSY partners of the ID infinite well are obtained. 
New potentials are thus produced that depend on two parameters (one is an integer and the 
other one is a real number). The complete energy spectra of the new Hamiltonians are obtained. 
The GCS are then constructed and their properties are exhibited. In Section 3, we generalise 
the preceding approach to 2D where the SUSY partners exhibit more freedom. Complete energy 
spectra are also produced. We generalise the definition of GCS to this context. We discuss as 
well another type of coherent states, the so-called generalised coherent states (GeCS). Because 
of the degeneracy of the energy spectrum, the usual definition of those states as the eigenstates 
of an annihilation operator has to be adjusted. We finish the paper with some conclusions and 
questions for future work. 


2. The ID infinite well and SUSY partners 

2.1. Description of the models 

Let us first set our notational convention concerning the infinite square well [1] to be used 
throughout this work. A particle of mass M is subject to a potential taken to be 


V(x) 


0 , 0 < X < IT 

oo, otherwise. 


(1) 


The stationary eigenstates and the discrete energies of this system are 


ipn(x) 



h 2 2 

E n = —— n , n = 1,2,... 
2M 


( 2 ) 


In the following, we will use dimensionless units, setting h = 1, M = 1/2, such that the 
Hamiltonian is H x = — + V(x). 

SUSY partners of the infinite well have been constructed m starting from usual intertwining 
relations involving the supercharges Q x , Q x . These can be defined, in particular, as differential 
operators of second order. Let us summarize the results. 

Starting from the Hamiltonian H x , a SUSY partner Hamiltonian H x = — fjf^ + V(x) is 
obtained from the relations 


fl,Q, = Q X H X , Q X H X = H x Ql 


( 3 ) 


For the infinite well, with V(x) = 0 in the domain x g]0,7t[, internal consistency constrains the 
supercharges to take the form 


and 


Q x = ^ + e + (y 2 ^) - v'(x)) 




( 4 ) 


( 5 ) 


where e is an arbitrary constant. The function r/(x) satisfies (in the so-called confluent case 

mr- 


2 r]{x)r)"{x) — ( rj'(x )) 2 — 4r/ 2 [x)r]' (x) + ij 4 (x) + 4 er] 2 (x) = 0 


(6) 


and the new potential is given as 


( 7 ) 


V(x) = 2 r/(x). 

Moreover, the products Q\Q X and Q X Q\ are respectively polynomials in H x and H x : 

QlQ x = (H x - e) 2 , Q x Ql = (H x -e) 2 . (8) 

The resolution of (|6| leads to admissible solutions for e = k 2 with k = 1,2,... A particular energy 
e = Ek within the original spectrum thus needs to be chosen for the supersymmetry. We get 

HH 


rj(x ; k, uj) = 


4 k sin 2 (kx) 


sin(2fex) + 2k(iru> — x) 
where uj is an arbitrary constant. The corresponding potentials are given as 

” 32k 2 sm(kx)[sin(kx)+k(nLJ—x) cos(fcrr)] c\ ^ / 

[s'm(2kx)-\-2k(7rLu—x )] 2 ’ ^ ^ 

oc, otherwise. 


V(x; k,cj) = 


(9) 


( 10 ) 


These potentials can be shown to be non singular if u> e] — oo,0[ U ]l,oo[. Two instances are 
illustrated on figure 1, one for each of the disconnected subsets of the parameter space of uj, i.e. 
uj = 2 and uj = —1. We thus observe that changing the sign of (uj — 1/2) reverses this potential 
about its vertical axis. Decreasing the magnitude of \uj — 1/2| accentuates the departure from a 
sinusoidal potential. The parameter k controls the amount of oscillations. 



Figure 1 - SUSY potential V(x\k= 10, uj) as a function of x, for (left) uj = 2 and (right) uj = — 1. 

Notice the striking symmetry under the Z 2 action w — > 1 — uj, x — > it — x. This symmetry 
is accidental in the sense that it is not a consequence of the SUSY algebra, nor the imposed 
intertwining relations <§• Still, it will have interesting effects on the coherent states built out 
of SUSY potential functions, as we shall discuss in section 2.3. 

The normalised SUSY eigenstates ipn( x ) are obtained from the intertwining relations ^ and 
the expressions Q: 

^n(x-, k,u) = (e - ^ n ) _1 Q x ^ n (x), (n^k). (11) 

More explicitly, we get 


^ n (x-,k,uj) 


2 sin(n.x)[sin(2fcx)(n 2 + k 2 ) + 2k(iruj — x)(n 2 


(n 2 — fc 2 )[sin(2foc) + 2 k( 


- k 2 )\ — 4nkcos(nx) sin 2 (kx) 
7 tuj — x)] 


(12) 























For n / k, they are physical states, i.e. they are normalisable and such that ip n (0]k,ui) = 
if n {i r; k,u) = 0, since 77 (0; k,u) = 77 ( 77 ; &, w) = 0. The corresponding energies are E n = n 2 as in 
the original case. 

For n = k, since e = E k = k 2 , this procedure does not yield if k (x) because Qx^kix) = 0. 
The completeness of the spectrum of H x has been investigated [H] and a single additional state 
was found. For the sake of completeness, we now explain how this missing state if k {x) could be 
obtained. It is, in fact, the simultaneous solution to the system of equations 

Qli>k(x) = 0, H x fi k (x) = e^fe(x), (13) 

that reduces to a first order differential equation on f>k{x)\ 


- rj(x;k,uj ) d + (V(x;k,uj) + ^(r/ 2 (x; k,u) - ))$ k (x) = 0. (14) 

From p3], we know that rj'(x) = r] 2 (x) + 2/3(x)r](x), with j3{x) solving the Riccati equation 
/3\x) + /3 2 (x) = —e. Equation (14) thus becomes 


dipkjx) 

dx 


= (r?(x; k, u) + /3(x- k))tp k (x) 


(15) 


and the normalised solution of energy e = E k = k 2 is given as 


Mx-Xu) = \/f Sin(fa) sin(2fa/+2*(J ) -x) - 


(16) 


It is normalisable and such that ip k (0, k,co) = k,co) = 0. With this additional state the 

spectrum of H x is thus complete. 


2.2. Gaussian Klauder coherent states for the infinite well 

As mentioned in the introduction, the gaussian Klauder coherent states (GCS) can be built 
for many different systems as a special superposition of energy eigenstates in order to get a 
reasonably well localised probability density distribution for a short period of time ( 9j. They 
have proven to be relevant for the study of the harmonic oscillator and the infinite well. For this 
last system, we recently m formalised the relation between those states and the generalised 
coherent states (GeCS), constructed as eigenstates of an annihilation operator of the system 
under consideration. We have shown in particular that the GCS can be chosen to reproduce 
approximately the GeCS in some specific area of their parameter space. In this section, we thus 
deal only with those GCS and summarise their properties. 

For real constants cj>o, no > 0 and <xo > 0, they are defined as the gaussian combination 

00 - in ~Zi — inct> ° 

^G(x,fin 0 ,a 0 ,(j) 0 ) = VC^(n o ,CT O ,c)>o)e _ * En Vn(30, (n 0 , 00 ,^ 0 ) = 6 7 = , (17) 

n= l V -^GWO; °o) 


where the normalisation factor is 

OO (n-nn ) 2 

77 G (n 0 ,cr 0 ) = 2<T o . (18) 

n=0 

The resolution of the identity is satisfied for theses states, as well as time stability and continuity 
in no and oq 0- 












Interestingly, not only © has a gaussian distribution in n, the sum can be carried out 
approximately to yield a gaussian wavefunction. Indeed, as demonstrated in 


V G(no,vo,<t>o-,x,t ) ~ 


(y/2n s ) 1 / 2 


exp 


(s ~ ^) 2 

4s 2 


+ iPx 


(19) 


(up to a x-independent phase factor), for x E [0, 7 r], t > 0 and with the following shorthands: 


Pt 


1 


X ^ 0+ l/2’ P n °’ * 2cr’ 4(r 2 + 1 2 ) ’ 


r = (40g) 


2 \ — 1 


Here stresses an approximate validity under the conditions no < 7 0 1, X s, ir—X S> s, 

and i«r. 

This result implies, as alluded to earlier, that the state is a well localised wave packet bouncing 
back and forth on the walls of the well, in a quasi-classical fashion. Its initial width is set by 
the parameter no, which also determines the time scale r of the Lorentzian decay of the packet. 
On the other hand, the velocity of the motion scales as the parameter no- Finally, the phase (j )o 
establishes the initial position of the wave packet. Here, of course, we naturally identify X ~ (x) 
and P ~ (p). 


Another remarkable consequence of (19) is that it explains why the GCS almost saturate 


Heisenberg uncertainty relation, when the conditions of validity are satisfied mi- This is readily 
understood by recalling the standard fact that the most general wave function minimizing Ax Ap 
is [5] 


T(x, t) = A(t)e 


~^~+i{p)x 


( 20 ) 


2.3. Gaussian Klauder coherent states for the SUSY partners 

The construction of the GCS can easily be attempted for the SUSY partner in just the same 
way it was performed in the original system: 


OO 


'F G (x,f;n o ,(Xo, 0 o) = V C^(n 0 , cr 0 , 0 o )e lEnt tf n (x) 


n= 1 


( 21 ) 


However, the algebraic construction reviewed in section 2.1 a priori does not enforce (21) to 


retain the well-understood quasi-classical behaviour summarised in (19). 


We notice on the other hand that the SUSY stationary states (11) with n 7 ^ k. can be 
regarded as perturbed versions of the initial eigenmodes ([ 2 ]), provided the overall sign is chosen 
as in (12). This is especially obvious at k 3> 1, i.e. by using a high energy mode for the SUSY 


mapping. 

The case n = k deserves here, like in section 2 . 1 , some special care. As clear from (16), the 
perturbative equivalence ~ i/; n (x'; k S> 1, w) no longer holds for this state. Thus, the mode 

ipk(x',k,(jo) corresponding to the energy used for the SUSY mapping is always fundamentally 
different from its original homologue ifk{x). 

The consequences of these remarks for the GCS now appear more distinctly. Whenever no is 
chosen far (i.e. out of a few standard deviations) from k. the GCS of the SUSY partner agree 
with the GCS of the original system. The analysis of section 2.2 then carries on, and a localised 
gaussian quasi-classical wave packet minimizing approximately Heisenberg relation is generated. 

If no and k are of the same order, on the other hand, some discrepant features are expected for 
the SUSY GCS. As a concrete example, figure 2 exhibits the time-evolution of the probability 
density distribution calculated from (21), in the coincident case no = k. Clearly, an effect 


of the supersymmetry is to generate superimposed wavelets on the principal gaussian wave 

















packet. As time evolves, moreover, some side modes appear next to the dominant packet. 
Numerical calculations show that these “harmonics” coherently follow the principal peak in its 
displacements inside the well as time evolves. The wave packet slowly spreads and bounces on 
the walls as in the original construction. 



(a) Initial distribution 



(b) Before 1st bounce 
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(c) After 1st bounce 


t« 0.012 




Figure 2 - Time-evolution of the probability density distribution inside the well for the SUSY coherent 
state ’F o{x,t;no = 100, ctq = 10, </>o = 7t/2). The SUSY parameters are: k = 100, u = 2. 


Let us now turn to the interesting question of the effect of the SUSY parameter u on the 
GCS. Our numerical calculations seem to suggest that the magnitude of |co — 1/21 weakly affect 
the behaviour of the states. Thus, the above discussion continues to apply. What is more subtle 
is the effect of changing the sign of (w — 1/2). In spite of its simple mirror effect on the potential, 
as discussed in section 2.1, it tremendously affects the GCS. As shown on figure 3, the Z 2 action 
suppresses the wavelets and side-modes, even for the coincident case no = k. This is surprising 
since the above discussion is still applicable, so we might still have expected some disagreement 
due to the exceptional state ipk(,x; k,u). At the level of the wave functions (12) and (16), the Z 2 
action acts as an alternating phase factor (except at energy Ek, where the alternating pattern 
is broken): i/; n -» — cos(mr)ilj n , ipk —> + cos(kn)'^k (provided we also map x —> tt — x). 

Note that although the GCS of the SUSY Hamiltonian can be chosen to be almost identical 
to the GCS of the original system, the potential functions remain very different. Indeed, the 
limit k 3> 1 emphatically does not lead to V(x;k,uj) — 0 in [0,7r], Thus, the construction of 
well behaved coherent states for this nontrivial system is a remarkable consequence of the SUSY 
approach. 





(a) Initial distribution 


(c) After 1st bounce 


(b) Before 1st bounce 

Figure 3 - Time-evolution of the probability density distribution inside the well for the SUSY coherent 
state 'S r G(x,t;no = 100, <jq = 10, d>o = tt/ 2). The SUSY parameters are: k = 100, uj = —1. 

As a concluding side comment, let us note that the classical dynamics in the SUSY 


system (10), which is generally expected for coherent state constructions, is not immediately 


transparent. A minimal criteria for classicality might thus be chosen to be a relatively slowly 























evolving and highly localized wave envelope. We have adhered to this point of view in the 
current work. 


3. The 2D infinite well and SUSY partners 

3.1. Description of the original model 

We consider again a particle of mass 1/2, now assumed to move in a square 2D box of size n. 
The Hamiltonian is given as 

H = H X + H y , (22) 

where H x and H y are the ID Hamiltonians of the infinite well in the directions x and y 
respectively. The corresponding normalised eigenstates and discrete energies are 

2 

'F n,m(x,y ) = if>n(x)il>m{.y) = - sin nx sin my (23) 

7r 


and 


E nm — ri fi TU . 


(24) 


where n,m = 1,2 ,... 

We see that this quadratic energy spectrum presents two types of degeneracies. The first type 
is a permutation degeneracy since the eigenstates y) and 4> min (x, y), n ^ m, are distinct 

with the same energy E ntfn . The second type is called accidental mm or arithmetic [16] . For 
example, we have E$ 5 = = 50 with 1 I' 5 i 5 (:r,y) and 4>ij(x,y) as distinct eigenfunctions 

and E\ t % = E^j = 65 with ’Fi } $(x,y) and d , 4,7 (x,y) as distinct eigenfunctions. The problem 
of identifying the number and type of these degeneracies for the infinite well has already been 
considered [nananE]. 

We present here some notations that will be useful in the discussion of coherent states of such 
a system: 


We rank the energies (24) in increasing order and write them as E v , with index v = 0,1,... 
We also introduce the shifted energies £ u = E u — Eq. 


• As alluded to above, more than one state may correspond to E u . We thus define fi v as 
the index for the summation of all the eigenstates associated to the degenerate energy E u 
and d v as the associated number of degeneracies. The index (jl v goes from 0 to d u — 1 . 
Although they are irrelevant for our purposes, let us point out that explicit formula for the 
existence and number of accidental degeneracies exist in the number theoretic literature. 
They involve the number of divisors of E u of the form 4k + 1 and 4k + 3, k G N. 

• The indices (n, m) corresponding to a state alternatively referred to with (v, y u ) are written 

• We rename the eigenstate )T7ll/ ^ (x, y), the state corresponding to the energy E v and 

degeneracy index /j, v , as (x,y). 

• To each energy E v . we can thus associate a single “cumulative” eigenstate constructed as a 
superposition of the states 4/ u, t J "'(x, y). It takes the form 


d v -1 

® u (x,y) = ^2 lv,nu i & v,tl “'{.x,y), (25) 

fj,u=0 

where the coefficients are arbitrary non zero complex numbers. Note that we have 
taken into account the permutation degeneracy in the superposition formula. 



3 


3,1 


P 3 = 
= arctan 



Figure 4 - Polar parametrisation of eigenstates 


• Finally, a polar parametrization {p u = y/Ef,, tan 9 U ^ V = is introduced. It provides 

an intuitive picture of the indices v and p v . As exemplified in figure 4, the energy index v is 
associated with the radius p u in (m, n)-plane, while the degeneracy index p v is associated 
with the angle 9 U ^ V . For definiteness, let us assume that p, u increases as n v increases, while 
m v decreases (for fixed v). Put differently, p u increases as the angle 0 V ^ V increases. 


3.2. Construction of SUSY partners 

2D generalisations of SUSY quantum mechanics have been proposed some time ago and, 
in particular, the approach of Ioffe and collaborators m has given some relevant results. 
Unfortunately, using supercharges that are essentially a linear combination of the ones used 
in ID does not lead to the new potentials we have already obtained in ID. An alternative 
approach is considered here. 

A SUSY partner Hamiltonian which is separable in x and y could be obtained from SUSY ID 
Hamiltonians as constructed in Section 2. Indeed, let us take Q = Q x Qy and Q t = QlQy where 
Qx, Qx and Q y , Qy satisfy the interwining ID relations §) in the corresponding variables. We 
thus easily get 

HQ^ = Q ] H, QH = HQ, (26) 


where H = H x + H y with H x = + V(x,k\,u\) and H y = + V(y, k 2 , w 2 ). The 

corresponding potentials are explicitly given in (10). 


d 2 

dy 2 


Eigenstates of the SUSY partner Hamiltonian H are thus easily constructed from 


n,m(x,y,ki, k 2 ,ui,u 2 ) = ((&? - n 2 )(*| - m 2 )) 1 Q + ^ n ,m(x,y) = ipn(x-, ki,ui )4>m{y, k 2 ,UJ 2 ), 


'P 

I V) I I u \ / tS / -*-/ / j- / / \\ X / \ ^ // W / *-/ / I • V \ / J. / J. / / \L/ / / 

(27) 

for n 7 ^ ki, m ^ k 2 . The corresponding energies are E n _ m = (24). In ID, we had to add one 
state (16) in order to get a complete spectrum. In 2D, we have to add more states. Indeed, we 
also get the following states: 

k 2 ,u 1,UJ 2 ) = v> fcl (x; fci, Wi)^ m (y; k 2 ,u 2 ), 

with energy E kum , 

vn,k 2 (x, y; ki, k 2 , Ui, UJ 2 ) = i} n {x\ ki, uJi)tp k2 (y; k 2 , u 2 ), 
with energy E n ^ 2 , and finally the eigenstate with energy E^ x ^ 2 is given as 

^kiM (x, y, ki, k 2 , Ui, u 2 ) = ip kl (x; h, iOi)^ k2 (y, k 2 , u 2 ). 


(28) 


(29) 


(30) 




















3.3. Coherent states 

Generalised and gaussian coherent states have been constructed for the usual 2D infinite well 
mu mi- The generalised coherent states (GeCS) are as usual defined as eigenstates of an 
annihilation operator of the system under consideration. But in the 2D case, due to the existence 
of a quadratic degenerate energy spectrum, we had to express the energies in increasing order 
and also to make a superposition of the different eigenstates with the same energy as given 
in (25). The construction of an annihilation operator and numerical calculations giving the 
behaviour of those states can be found in m- They have been shown to be closely related to 
GCS for some values of the coherent states parameters. 

We will slightly adjust these constructions in this work. 


3.3.1. Gaussian coherent states The definition of GCS is proposed as a direct generalisation 
of the ID case. Indeed, we take 

^ 2 ( P(x,y,t-,n 0 ,m 0 ,a no ,a mo ,(j) no ,(j) mo ) = T G (.T,f;n 0 i @no i 0no )^G(y,t-,m 0 ,cr mo ,(j)mo), (31) 

with Tq = ©■ Explicitly, we get 


OO OO 

^{x,y,t;no,m 0 ,a no ,(Tm 0 ,(t>n 0 ,(l)m 0 ) = C n,m e ~ iEn ’ mt ^m,n( x ^y)^ ( 32 ) 

n =1 m =1 


where 


Cn,m — ^'n,m{ n 0t m 0i (T no , <X mo , (f> no , — C 71 (no, a no , (j) no )C m (m 0 , a mo , (j) mo ). (33) 

This definition pays off since the analysis of El reviewed in section 2.2 continues to apply 
here. As in the ID case, we get a well localised state and quasi-classical behaviour. Approximate 
saturation of the Heisenberg uncertainty bound is also achieved, as in the ID setting. 



Figure 5 - Probability density distribution and time evolution of the position expectation value of the 
SUSY GCS in 2D until t = 0.012. State parameters are chosen as: no = mg = 100, a no = a mo = 10, 
(f>no = 7 t/2 , and <j> no = 0. SUSY parameters are: k\ = 100, &2 = 50, wi = —1, and lo 2 = 2. 






Such a definition of GCS is easily extended to the SUSY case where the eigenstates *1 ’ m ,n(x, y ) 

i(x,y). Again, the factorization of the state yields 
As an example, 


in (32) are replaced by the eigenstates v b r 


a straightforward generalisation of the ID results. As an example, figure 5 shows a trace of 
the time evolution of the position expectation value for the SUSY GCS in 2D. A neat localized 
packet bouncing on the walls is obtained once again. 


3.3.2. Generalised coherent states The GeCs are eigenstates of an annihilation operator of the 
2D infinite well and they have been defined as m 


^Ge(x,y,t-z) 


E 


—i£vt 


V N Ge(z) ^ y/p{y) 


$ u (x,y), 


(34) 


with the normalisation factor 


^Ge(*) = £ 


°° 'z\ 2v 


PM 


and 


P( v ) 


1 , if v = 0 

nr=i^, if^^o. 


(35) 


These states depend on a continuous complex parameter z as in the ID case. With respect to 
an alternative approach |12| . more freedom is given to the states & u (x,y) being a superposition 
of states with same energy £ v . 


4. Conclusion and future work 

SUSY partners of the infinite well have been constructed and have shown to satisfy relevant 
properties compared with the ones of the original quantum system. In particular, for the 
ID system, we have constructed a set of coherent states (GCS) which depends on real and 
discrete parameters. The relation between those states and the generalised ones (GeCS) has 
been formally given m and we have shown the behaviour of the GCS case. 

In the 2D case, we have extended the construction of SUSY partners to get similar potentials 
as in ID. More parameters are involved in this context and the existence of degeneracies in the 
energy spectrum has lead us to adjust the definition of coherent sates with respect with the ones 
used in the ID case. The GCS are constructed in order to have a good behaviour with respect 
to localisation in the usual case as well as in the SUSY case. We have also constructed GeCS in 
the 2D setting. 

It remains to make a link between both types of states. The use of the polar parametrisation 
is clearly a way to solve the problem. Indeed, the GCS may be written as 


^c(x,y,t) = 


°° (Ph'—Pp ) 2 




E' 


0 iSjy t lTi G, IS 


^ ,v {x,y). 


(36) 


The gaussian terms depending on the radial variables could thus be related to the factor 
z u {p{v))~ 2 in (34) as it was the case in ID [11] . The gaussian terms involving the angular 
variables are related to the coefficients of the superposition of eigenstates with same energy. 
This way, the states ^ G,v {x,y) will be related to $> v (x,y) in (34). As mentioned earlier, a 
particular approach of this question in 2D mm has shown that GeCS are in fact a good 
approximation of GCS. In a future work, we hope to formally solve this problem. 
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